Keratectomy data determining device and keratectomy data determining program

ABSTRACT

The purpose of the present invention is to provide a corneal-ablation-data determining apparatus and a corneal-ablation-data determining program capable of controlling corneal ablation amount distribution and expanding ablation applicability. In the present invention, the corneal-ablation-data determining apparatus for use in keratorefractive surgery for ablating a corneal surface to correct a refractive error has input means for inputting measurement data on any one of wavefront aberration distribution and refractive power distribution of a patient&#39;s eye, decentering means for decentering a reference axis being a calculation reference for a corneal ablation amount, and calculation means for obtaining the corneal ablation amount based on the measurement data and data on eccentricity of the reference axis.

TECHNICAL FIELD

[0001] The present invention relates to a corneal-ablation-data determining apparatus and a corneal-ablation-data determining program for use in keratorefractive surgery for ablating a corneal surface of a patient's eye (an eye to be operated on) to correct a refractive error.

BACKGROUND ART

[0002] There is known keratorefractive surgery for ablating a corneal surface with a laser beam in an ultraviolet wavelength range and changing its shape to correct a refractive error of a patient's eye. In a case where even aberration of the eye is to be corrected by the surgery, mainly two methods are used for calculation of a corneal ablation amount.

[0003] The first one is a method of obtaining the corneal ablation amount based on measurement data on wavefront aberration in a wide area of the patient's eye. In this method, the corneal ablation amount may be obtained using Expression. 1 being an approximate expression which is conventionally known as a processing technique for an optical component. $\begin{matrix} {P \approx \frac{\delta \quad W}{n - 1}} & \text{[Expression~~1]} \end{matrix}$

[0004] In Expression 1, δW is wavefront aberration, n is a refractive index of an optical component to be corrected, P is a correction amount for a surface.

[0005] The second one is, as described in U.S. Pat. No. 6,033,075 (Japanese Patent Application Unexamined Publication No. Hei11-342152) and U.S. Pat. No. 6,086,204, a method of obtaining the corneal ablation amount using measurement data on objective refractive power distribution or that on wavefront aberration distribution in the wide area of the patient's eye, and data on a pre-operative corneal shape.

[0006] The corneal ablation amount is obtained by these conventional calculation methods as an ablation amount where the patient's eye is emmetropia with no aberration in a measurement optical axis direction established at the time of the measurement of the objective refractive power distribution or the wavefront aberration distribution (or a non-aberration eye which has particular refractive power in the measurement optical axis direction). However, a human eye is an eccentric optical system comprised of refractive elements such as the cornea and a crystalline lens, and the corneal shape is generally asymmetrical with respect to a visual axis. Therefore, in the conventional calculation methods for the corneal ablation amount, there is a problem that the ablation amount becomes eccentric.

[0007] In addition, in the case of additional ablation performed for correcting the patient's eye after the ablation through axial-deviated irradiation of the laser beam, there is a case where a proper corneal ablation amount cannot be obtained with the conventional calculation method for the corneal ablation amount.

[0008] The present invention has been made in view of the above circumstances and has an object to overcome the above problems and to provide a corneal-ablation-data determining apparatus and a corneal-ablation-data determining program capable of controlling corneal ablation amount distribution and expanding ablation applicability.

DISCLOSURE OF THE INVENTION

[0009] To achieve the objects and in accordance with the purpose of the present invention, as embodied and broadly described herein, the present invention has the following features.

[0010] 1) A corneal-ablation-data determining apparatus for use in keratorefractive surgery for ablating a corneal surface to correct a refractive error is characterized as having input means for inputting measurement data on any one of wavefront aberration distribution and refractive power distribution of a patient's eye, decentering means for decentering a reference axis being a calculation reference for a corneal ablation amount, and calculation means for obtaining the corneal ablation amount based on the measurement data and data on eccentricity of the reference axis.

[0011] 2) In the corneal-ablation-data determining apparatus according to 1), the decentering means includes means for decentering the reference axis with respect to a measurement optical axis used for measurement of any one of the wavefront aberration distribution and the refractive power distribution, and means for inputting the data on the eccentricity.

[0012] 3) In the corneal-ablation-data determining apparatus according to 1), the input means inputs the measurement data on the refractive power distribution, and the calculation means obtains a post-operative corneal shape based on the measurement data on the refractive power distribution using a vector expression under the Snell's law so that an optical path inside the eye becomes the same before and after the surgery.

[0013] 4) In the corneal-ablation-data determining apparatus according to 1) characterized as further having means for inputting data on a pre-operative corneal shape, the calculation means obtains a post-operative corneal shape based on the data on the pre-operative corneal shape, the measurement data, and the data on the eccentricity of the reference axis.

[0014] 5) In the corneal-ablation-data determining apparatus according to 1), the calculation means obtains a direction vector of incident light into a post-operative cornea when the reference axis is decentered, and obtains a condition for a post-operative corneal shape as a normal vector of the cornea based on the direction vector.

[0015] 6) A corneal-ablation-data determining program for calculating a corneal ablation amount for use in keratorefractive surgery for ablating a corneal surface to correct a refractive error, the program is characterized as making a computer to function as means for inputting measurement data on any one of wavefront aberration distribution and refractive power distribution of a patient's eye, means for inputting data on eccentricity for decentering a reference axis being a calculation reference for a corneal ablation amount with respect to a measurement optical axis used for the measurement, and calculation means for obtaining the corneal ablation amount based on the measurement data and the data on the eccentricity of the reference axis.

[0016] 7) In the corneal-ablation-data determining program according to 6), the calculation means obtains a direction vector of incident light into a post-operative cornea when the reference axis is decentered, and obtains a condition for a post-operative corneal shape as a normal vector of the cornea based on the direction vector.

BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments of the present invention and, together with the description, serve to explain the objects, advantages and principles of the invention. In the drawings,

[0018]FIG. 1 is a view showing a schematic configuration of an optical system in a corneal-ablation-amount determining apparatus consistent with the present invention;

[0019]FIG. 2 is a view showing an arrangement of photodetectors in an eye refractive power measurement optical system;

[0020]FIG. 3 is a view showing a block diagram of a control system in the corneal-ablation-amount determining apparatus;

[0021]FIG. 4 is a view illustrating configuration of a coordinate system when a post-operative reference axis is decentered;

[0022]FIG. 5 is a view illustrating an example of a display of a simulation map at the time of calculation of a corneal ablation amount; and

[0023]FIG. 6 is a view describing a case where additional ablation is performed on a patient's eye after ablation through axial-deviated irradiation.

BEST MODE FOR CARRYING OUT THE INVENTION

[0024] A detailed description of one preferred embodiment of the present invention will now be given referring to the accompanying drawings. FIG. 1 is a view showing a schematic configuration of an optical system in a corneal-ablation-amount determining apparatus consistent with the present invention. The optical system is roughly divided into an eye refractive power measurement optical system, a fixation target optical system and a corneal curvature measurement optical system.

[0025] (Eye Refractive Power Measurement Optical System)

[0026] Reference numeral 100 indicates the eye refractive power measurement optical system consisting of a slit light projection optical system 1 and a slit image detection optical system 10. Near infrared light emitted from a light source 2 of the projection optical system 1 is reflected by a mirror 3 to illuminate a slit opening 4 a in a rotation sector 4. A motor 5 rotates the sector 4. Slit light scanned by the rotation of the sector 4 passes through a projection lens 6 and a limit diaphragm 7, and is reflected by a beam splitter 8. Then, the slit light is transmitted through a beam splitter 9 which makes an optical axis of the fixation target optical system and an observation optical system (described later) coaxial, and converges in the vicinity of a cornea Ec of a patient's eye (an eye to be examined, an eye to be operated on) E to be projected onto a fundus Ef. Besides, the light source 2 is arranged at a position conjugate with the vicinity of the cornea Ec with respect to the lens 6.

[0027] The detection optical system 10 is provided with a photo-receiving lens 11 and a mirror 12 disposed on a principal optical axis L1, and a diaphragm 13 and a photo-receiving unit 14 disposed on an optical axis L3, which is formed by reflection of the mirror 12. The diaphragm 13 is arranged at a back focal point of the lens 11 via the mirror 12. (That is to say, it is arranged at a position conjugate with an emmetropic fundus.) As shown in FIG. 2, the photo-receiving unit 14 includes, on its photo-receiving surface, eight photodetectors 15 a-15 h arranged at positions approximately conjugate with the cornea Ec with respect to the lens 11. Among them, the photodetectors 15 a-15 f are positioned on a line passing through a center of the photo-receiving surface (the optical axis L3), and the photodetectors 15 a and 15 b, 15 c and 15 d, 15 e and 15 f, respectively, are arranged symmetrically with respect to the center of the photo-receiving surface. As for these three pairs of the photodetectors, distance in the arrangement is set so as to detect refractive power corresponding to each position in a meridian direction on the cornea Ec. (In FIG. 2, it is shown as an equivalent size on the cornea.) On the other hand, the photodetectors 15 g and 15 h are arranged symmetrically with respect to the optical axis L3 on a line perpendicular to the line on which the photodetectors 15 a-15 f are arranged.

[0028] In the measurement optical system 100 having a constitution described above, a rotation mechanism 21 consisting of a motor 20, a gear and the like rotates synchronously the light source 2 to the motor 5 in the projection optical system 1 about the optical axis L2, and the photo-receiving unit 14 about the optical axis L3. The preferred embodiment is configured so that, when the slit light is scanned on a fundus of a hyperopic eye or a myopic eye without astigmatism, the arrangement direction of the photodetectors 15 a-15 f is perpendicular to a longitudinal direction of the slit light (image) photo-received on the photo-receiving unit 14.

[0029] (Fixation Target Optical System)

[0030] Reference numeral 30 is the fixation target optical system, 31 is a visible light source, 32 is a fixation target, and 33 is a projection lens. The lens 33 moves in an optical axis direction to fog the eye E. A beam splitter 34 makes the optical axis of the observation optical system coaxial. The light source 31 illuminates the fixation target 32. Light from the fixation target 32 passes through the lens 33 and the beam splitter 34, and is reflected by the beam splitter 9 to head for the eye E. Thereby, the eye E may fixate the fixation target 32.

[0031] (Corneal Radius of Curvature Measurement Optical System)

[0032] The corneal radius of curvature measurement optical system consists of a target projection optical system 25 for measuring a corneal radius of curvature and a target detection optical system 35 for measuring a corneal radius of curvature. The projection optical system 25 has the following constitution. Reference numeral 26 is a conic placido plate which has an opening in a center part. Formed on the placido plate 26 is a ring-pattern, which has a number of light transmission parts and light shielding parts concentrically arranged while having the optical axis L1 as the center. Reference numeral 27 is a plurality of illumination light sources such as an LED. Illumination light emitted from the light sources 27 is reflected by a reflection plate 28 to illuminate the placido plate 26 from behind almost uniformly. The light of the ring-pattern transmitted through the transmission parts of the placido plate 26 is projected onto the cornea Ec to form a ring-pattern (placido ring) image thereon.

[0033] The detection optical system 35 is provided with the beam splitter 9, the beam splitter 34, a photographing lens 37 and a CCD camera 38. The light of the ring-pattern image formed on the cornea Ec is reflected by the beam splitters 9 and 34, and enters (is received by) an image-pickup element of the camera 38 by the lens 37. In addition, the detection optical system 35 doubles as the observation optical system, and the light of an anterior-segment image of the eye E illuminated by an anterior-segment illumination light source unillustrated enters (is received by) the image-pickup element of the camera 38. A TV monitor 39 displays the photographed anterior-segment image and the ring-pattern image.

[0034]FIG. 3 is a block diagram of a control system. A measurement system control unit 50 is connected with the photodetectors 15 a-15 h, the camera 38, the monitor 39, a measurement switch 54, a memory 55, an FDD (floppy disk drive) 56, and the like to control each element. The monitor 39 is formed of a touch panel, and doubles as input means. The control unit 50 is provided with a refractive power calculation unit 52 and a corneal shape calculation unit 53. Measurement data is recorded on a floppy disk 59 by the FDD 56, and inputted to a computer 60. The computer 60 is provided with a control unit 61, a color display 62, a keyboard 63, a mouse 64, an FDD 66 and the like. The control unit 61 includes a calculation unit 61 a and a memory 61 b which stores a control program.

[0035] Next, measurement of eye refractive power distribution and a corneal radius of curvature (a corneal shape) will be described hereinafter. A touch key on the monitor 39 is operated to place the apparatus in an eye refractive power measurement mode. An examiner performs alignment while observing the image of the anterior-segment of the eye E on the monitor 39. As for the alignment, a target for alignment unillustrated is projected onto the cornea Ec from the front of the eye E so that its corneal reflex and a reticle may have a predetermined positional relationship. In the apparatus in the preferred embodiment, the measurement optical axis (L1) is aligned with an axis connecting the corneal reflex and the fixation target. After the completion of the alignment, the measurement switch 54 is depressed and a trigger signal is generated to perform the measurement.

[0036] The eye refractive power calculation unit 52 obtains objective eye refractive power distribution based on phase difference of an output signal from each of the photodetectors in the photo-receiving unit 14. First, preliminary measurement is performed in the same manner as refractive power measurement in a conventional phase difference method, and the lens 33 is moved based on its result to fog the eye E. After that, a corneal vertex (a corneal center) in a meridian direction on which the photodetectors 15 a-15 f are placed is obtained based on the output signals from the photodetectors 15 g and 15 h which vary in accordance with movement of the slit light (image) on the photo-receiving part 14. Next, the refractive power in a corneal portion corresponding to each of the photodetectors is obtained based on the phase difference of the output signal from each of the photodetectors 15 a-15 f with respect to the corneal vertex (corneal center). Then, while making the projection optical system 1 and the photo-receiving unit 14 to rotate 180 degrees at a step with a predetermined angle (1 degree), the refractive power is calculated for every meridian at each step of angles to obtain the distribution of the refractive power which varies in the meridian direction. (For more details, please refer to U.S. Pat. No. 5,907,388 (Japanese Patent Application Unexamined Publication No. Hei10-108837).) Here, an eye refractive power value is obtained while taking the corneal vertex (corneal center) as reference; however, it may also be obtained while taking a spectacle-lens wearing position as the reference. Data on the objective eye refractive power obtained is stored in the memory 55.

[0037] At the time of the measurement of the corneal radius of curvature, the touch key on the monitor 39 is operated to place the apparatus in a corneal radius of curvature measurement mode. The alignment is performed in the same manner as in the eye refractive power measurement mode, and the measurement switch 54 is depressed to perform measurement. The corneal shape calculation unit 53 processes the image photographed by the camera 38 to detect an edge of the ring-pattern image. Then, each edge position with respect to the corneal vertex (corneal center) is obtained for every step of the predetermined angle (1 degree) to calculate the corneal radius of curvature. (For more details, please refer to U.S. Pat. No. 5,500,697 (Japanese Patent Application Unexamined Publication Hei7-124113) and others.) Data on the corneal radius of curvature obtained is stored in the memory 55.

[0038] When the measurement data on the objective eye refractive power and that on the corneal radius of curvature are obtained as provided above, they are inputted to the computer 60 via the floppy disk 59, and a corneal ablation amount is calculated by operating the keyboard 63, the mouse 64 and the like provided to the computer 60. Hereinafter, a calculation method of the corneal ablation amount will be described.

[0039] <Condition of Post-Operative Corneal Shape Derived from the Snell's Law>

[0040] In a case where even the aberration is corrected by the corneal ablation, it is necessary that an optical path inside the eye becomes the same before and after the ablation. Therefore, inclination of each point on a corneal surface is varied accordingly to absorb variation in an angle of incidence into the cornea. When this is expressed using Formula 1 being a vector formula under the Snell's law, Formulae 2 and 3 are obtained. $\begin{matrix} {{\left\lbrack {{Expression}\quad 2} \right\rbrack n^{\prime}\overset{->}{Q^{\prime}} \times \overset{->}{N}} = {n\overset{->}{Q} \times \overset{->}{N}}} & {{Formula}\quad 1} \\ {{\left\lbrack {{Expression}\quad 3} \right\rbrack {{n_{c}{\overset{->}{Q}}_{c} \times {\overset{->}{N}}_{pre}} = {{\overset{->}{Q}}_{pre} \times {\overset{->}{N}}_{pre}}}}{and}} & {{Formula}\quad 2} \\ {{n_{c}{\overset{->}{Q}}_{c} \times {\overset{->}{N}}_{post}} = {{\overset{->}{Q}}_{post} \times {\overset{->}{N}}_{post}}} & {{Formula}\quad 3} \end{matrix}$

[0041] n_(c): a refractive index of the cornea

[0042] {right arrow over (Q)}_(c): a unit direction vector of refractive light in the cornea

[0043] {right arrow over (Q)}_(pre): a unit direction vector of incident light into the cornea before the correction

[0044] {right arrow over (Q)}_(post): a unit direction vector of the incident light into the cornea after the correction

[0045] {right arrow over (N)}_(pre): a unit normal vector of the cornea before the correction

[0046] {right arrow over (N)}_(post): a unit normal vector of the cornea after the correction

[0047] The vector N_(pre) may be known by the corneal shape measurement. In addition, the vector Q_(pre) may be obtained as a unit normal vector of wavefront aberration into which the measurement result of the above-mentioned refractive power distribution is once converted. In general, abeam emitted from a central fovea going out of the eye does not necessarily cross the measurement optical axis, because a meridian direction component is obtained in the measurement of the refractive power distribution. In order to calculate the post-operative corneal shape to obtain the vector Q_(post) which is designated by using the vector N_(pre) and the vector Q_(pre), the following process is to be taken.

[0048] First, a vector n_(c)Q_(c) is obtained using Formula 4 which is equivalent to Formula 1. $\begin{matrix} {\left\lbrack {{Expression}\quad 4} \right\rbrack {{n^{\prime}\overset{->}{Q^{\prime}}} = {{n\overset{->}{Q}} + {\left\{ {{\pm \sqrt{n^{\prime 2} - n^{2} + {n^{2}\left( {\overset{->}{Q} \cdot \overset{->}{N}} \right)}^{2}}} - {\overset{->}{Q} \cdot \overset{->}{N}}} \right\} \overset{->}{N\quad}\left\{ \begin{matrix} {{{+ \cdots}\quad {\overset{->}{Q} \cdot \overset{->}{N}}} \geq 0} \\ {{{- \cdots}\quad {\overset{->}{Q} \cdot \overset{->}{N}}} < 0} \end{matrix} \right.}}}} & {{Formula}\quad 4} \end{matrix}$

[0049] In the case of Formula 2, Formula 5 is given. $\begin{matrix} {\left\lbrack {{Expression}\quad 5} \right\rbrack {n_{c}{\overset{->}{Q}}_{c}} = {{\overset{->}{Q}}_{pre} + {\left\{ {{\pm \sqrt{n_{c}^{2} - 1 + \left( {{\overset{->}{Q}}_{pre} \cdot {\overset{->}{N}}_{pre}} \right)^{2}}} - {{\overset{->}{Q}}_{pre} \cdot {\overset{->}{N}}_{pre}}} \right\} {\overset{->}{N}}_{pre}\left\{ \begin{matrix} {{{+ \cdots}\quad {{\overset{->}{Q}}_{pre} \cdot {\overset{->}{N}}_{pre}}} \geq 0} \\ {{{- \cdots}\quad {{\overset{->}{Q}}_{pre} \cdot {\overset{->}{N}}_{pre}}} < 0} \end{matrix} \right.}}} & {{Formula}\quad 5} \end{matrix}$

[0050] Thereby, the vector nQ_(c) is determined.

[0051] Next, the vector N_(post) is obtained by Formula 3. Since the vector nQ_(c) is already known and the vector Q_(post) is a value to be designated, Formula 3 may be solved for the unknown vector N_(post). Formula 3 is transformed to Formula 6. $\begin{matrix} {\left\lbrack {{Expression}\quad 6} \right\rbrack {{\left( {{n_{c}{\overset{->}{Q}}_{c}} - {\overset{->}{Q}}_{post}} \right) \times {\overset{->}{N}}_{post}} = 0}} & {{Formula}\quad 6} \end{matrix}$

[0052] As a vector product of two vectors (vectors that are not zero) is zero, they are parallel to each other. Also, the vector N_(post) is a unit vector. Therefore, the vector N_(post) may be expressed as Formula 7. $\begin{matrix} {\left\lbrack {{Expression}\quad 7} \right\rbrack {{\overset{->}{N}}_{post} = \frac{{n_{c}{\overset{->}{Q}}_{c}} - {\overset{->}{Q}}_{post}}{{{n_{c}{\overset{->}{Q}}_{c}} - {\overset{->}{Q}}_{post}}}}} & {{Formula}\quad 7} \end{matrix}$

[0053] If this is performed for each point using Formulae 5 and 7, the condition of the post-operative corneal shape may be obtained as the normal vector at each point.

[0054] Incidentally, in the present embodiment, the post-operative corneal shape is calculated so as to correct the aberration of the eye based on the measurement result of the refractive power distribution. Since the method for formulating the calculation condition for the post-operative corneal shape using the vector formula under the Snell's law as above does not perform paraxial calculation, more precise post-operative corneal shape may be obtained.

[0055] <Method for Obtaining a Unit Direction Vector of the Incident Light Into the Post-Operative Cornea (Q_(post)) When Eccentricity is Applied to the Post-Operative Reference Axis>

[0056] In order to use the method for obtaining the vector N_(post) as above, the vector Q_(post) is expressed as an expression. The measurement optical axis of the eye refractive power measurement (and the corneal shape measurement) is set as reference coordinates, and post-operative reference axis is assumed to be applied with eccentricity in shift and tilt from the reference coordinates.

[0057]FIG. 4 is a view illustrating configuration of a coordinate system when the post-operative reference axis is decentered. In FIG. 4, an axis OZ is the measurement optical axis, and an axis PO′ is the post-operative reference axis. As for the eccentricity of the post-operative reference axis PO′, the shift is designated by δx and δy, and the tilt is designated by φ and φ. In addition, assuming that the post-operative refractive power (Refractive Error) is targeted to become D[diopter] (a spherical value), Formula 8 is presented. $\begin{matrix} {\left\lbrack {{Expression}\quad 8} \right\rbrack {\overset{\rightarrow}{O^{\prime}P} = {{- \frac{1000}{D}}\begin{pmatrix} {\sin \quad {\varphi cos}\quad \phi} \\ {\sin \quad {\varphi sin}\quad \phi} \\ {{- \cos}\quad \varphi} \end{pmatrix}}}} & {{Formula}\quad 8} \end{matrix}$

[0058] Here, if the vectors Q_(pre) and Q_(post) for a point C on the cornea is considered, the vector Q_(post) being a unit direction vector parallel to a line PC is firstly expressed as Formula 9. $\begin{matrix} {\left\lbrack {{Expression}\quad 9} \right\rbrack \begin{matrix} {{\overset{->}{Q}}_{post} = \frac{\overset{\rightarrow}{PC}}{\overset{\rightarrow}{PC}}} \\ {= \frac{\begin{pmatrix} {x - {\delta \quad x} + {\left( {1000\quad \sin \quad {\varphi cos}\quad \phi} \right)/D}} \\ {y - {\delta \quad y} + {\left( {1000\quad \sin \quad {\varphi sin}\quad \phi} \right)/D}} \\ {{f_{c}\left( {x,y} \right)} - {\left( {1000\quad \cos \quad \varphi} \right)/D}} \end{pmatrix}}{\sqrt{\begin{matrix} {\left\{ {x - {\delta \quad x} + {\left( {1000\quad \sin \quad {\varphi cos}\quad \phi} \right)/D}} \right\}^{2} +} \\ \begin{matrix} {\left\{ {y - {\delta \quad y} + {\left( {1000\quad \sin \quad {\varphi sin}\quad \phi} \right)/D}} \right\}^{2} +} \\ \left\{ {{f_{c}\left( {x,y} \right)} - {\left( {1000\quad \cos \quad \varphi} \right)/D}} \right\}^{2} \end{matrix} \end{matrix}}}} \\ {= {\frac{\frac{D}{D}}{\sqrt{\begin{matrix} {\left\{ {\frac{\left( {x - {\delta \quad x}} \right)D}{1000} + {\sin \quad {\varphi cos}\quad \phi}} \right\}^{2} +} \\ {\left\{ {\frac{\left( {y - {\delta \quad y}} \right)D}{1000} + {\sin \quad {\varphi sin}\quad \phi}} \right\}^{2} +} \\ \left\{ {\frac{{f_{c}\left( {x,y} \right)}D}{1000} - {\cos \quad \varphi}} \right\}^{2} \end{matrix}}}\begin{pmatrix} {\frac{\left( {x - {\delta \quad x}} \right)D}{1000} + {\sin \quad {\varphi cos\phi}}} \\ {\frac{\left( {y - {\delta \quad y}} \right)D}{1000} + {\sin \quad {\varphi sin}\quad \phi}} \\ {\frac{{f_{c}\left( {x,y} \right)}D}{1000} - {\cos \quad \varphi}} \end{pmatrix}}} \end{matrix}} & {{Formula}\quad 9} \end{matrix}$

[0059] <Method for Obtaining the Post-Operative Corneal Shape>

[0060] Next, a method will be described which is for obtaining the post-operative corneal shape based on the normal vector of the post-operative cornea obtained by applying Formulae 5 and 7 to all data points.

[0061] First, the normal vector N when the corneal shape is given in an x-y coordinates form (see FIG. 4) is formulated. In this case, a position vector r at each position is expressed as Formula 10.

[0062] [Expression 10]

{right arrow over (r)}=x{right arrow over (i)}+y{right arrow over (j)}+f _(c)(x,y){right arrow over (k)}  Formula 10

[0063] By partially differentiating the above formula with respect to x and y, respectively, Formula 11 is given. $\begin{matrix} {\left\lbrack {{Expression}\quad 11} \right\rbrack {\frac{\partial\overset{->}{r}}{\partial x} = {\overset{->}{i} + {\frac{\partial{f_{c}\left( {x,y} \right)}}{\partial x}\overset{->}{k}}}},{\frac{\partial\overset{->}{r}}{\partial y} = {\overset{->}{j} + {\frac{\partial{f_{c}\left( {x,y} \right)}}{\partial y}\overset{->}{k}}}}} & {{Formula}\quad 11} \end{matrix}$

[0064] Accordingly, using Formula 12, the normal vector N is obtained by Formula 13. $\begin{matrix} {\left\lbrack {{Expression}\quad 12} \right\rbrack \begin{matrix} {{\frac{\partial\overset{->}{r}}{\partial x} \times \frac{\partial\overset{->}{r}}{\partial y}} = {\begin{matrix} \overset{->}{i} & \overset{->}{j} & \overset{->}{k} \\ 1 & 0 & \frac{\partial{f_{c}\left( {x,y} \right)}}{\partial x} \\ 0 & 1 & \frac{\partial{f_{c}\left( {x,y} \right)}}{\partial y} \end{matrix}}} \\ {= {{{- \frac{\partial{f_{c}\left( {x,y} \right)}}{\partial x}}\overset{->}{i}} - {\frac{\partial{f_{c}\left( {x,y} \right)}}{\partial y}\overset{->}{j}} + \overset{->}{k}}} \end{matrix}} & {{Formula}\quad 12} \end{matrix}$

$\begin{matrix} {\left\lbrack {{Expression}\quad 13} \right\rbrack \begin{matrix} {{\overset{->}{N} = {\frac{\partial\overset{->}{r}}{\partial x} \times {\frac{\partial\overset{->}{r}}{\partial y}/{{\frac{\partial\overset{->}{r}}{\partial x} \times \frac{\partial\overset{->}{r}}{\partial y}}}}}}\quad} \\ {= {\frac{1}{\sqrt{\left( \frac{\partial{f_{c}\left( {x,y} \right)}}{\partial x} \right)^{2} + \left( \frac{\partial{f_{c}\left( {x,y} \right)}}{\partial y} \right)^{2} + 1}}\begin{pmatrix} {- \frac{\partial{f_{c}\left( {x,y} \right)}}{\partial x}} \\ {- \frac{\partial{f_{c}\left( {x,y} \right)}}{\partial y}} \\ 1 \end{pmatrix}}} \end{matrix}} & {{Formula}\quad 13} \end{matrix}$

[0065] When the corneal shape is given as a polar coordinate system (ρcosθ, ρsinθ, fp(ρ,θ))(see FIG. 4), then Formula 14 is presented.

[0066] [Expression 14]

{right arrow over (r)}=ρ cos θ{right arrow over (i)}+ρ sin θ{right arrow over (j)}+f _(p)(ρ,θ){right arrow over (k)}  Formula 14

[0067] {right arrow over (r)}: a position vector

[0068] {right arrow over (i)} a unit vector in the x direction

[0069] {right arrow over (j)} a unit vector in the y direction

[0070] {right arrow over (k)} a unit vector in the z direction

[0071] Accordingly, Formula 15 is obtained. $\begin{matrix} {\left\lbrack {{Expression}\quad 15} \right\rbrack {{\frac{\partial\overset{}{r}}{\partial\rho} = {{\cos \quad \theta \quad \overset{}{i}} + {\sin \quad \theta \quad \overset{}{j}} + {\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}\overset{}{k}}}},{\frac{\partial\overset{}{r}}{\partial\theta} = {{{- {\rho sin}}\quad \theta \quad \overset{}{i}} + {{\rho cos}\quad \theta \quad \overset{}{j}} + {\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta}\overset{}{k}}}}}} & {{Formula}\quad 15} \end{matrix}$

[0072] From Formula 15, Formula 16 is given. $\begin{matrix} {\left\lbrack {{Expression}\quad 16} \right\rbrack {{\frac{\partial\overset{}{r}}{\partial\rho} \times \frac{\partial\overset{}{r}}{\partial\theta}} = {{\begin{matrix} \overset{}{i} & \overset{}{j} & \overset{}{k} \\ {\cos \quad \theta} & {\sin \quad \theta} & \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho} \\ {{- {\rho sin}}\quad \theta} & {\rho \quad \cos \quad \theta} & \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta} \end{matrix}} = {{\left( {{\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta}\sin \quad \theta} - {\rho \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}\cos \quad \theta}} \right)\overset{}{i}} - {\left( {{\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta}\cos \quad \theta} + {\rho \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}\sin \quad \theta}} \right)\overset{}{j}} + {\rho \overset{}{k}}}}}} & {{Formula}\quad 16} \end{matrix}$

[0073] Therefore, the normal vector N is expressed as Formula 17. $\begin{matrix} {\left\lbrack {{Expression}\quad 17} \right\rbrack {\overset{}{N} = {{\frac{\partial\overset{}{r}}{\partial\rho} \times {\frac{\partial\overset{}{r}}{\partial\theta}/{{\frac{\partial\overset{}{r}}{\partial\rho} \times \frac{\partial\overset{}{r}}{\partial\theta}}}}} = {\cdots = {\frac{1}{\sqrt{\left( {\rho \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}} \right)^{2} + \left( \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta} \right)^{2} + \rho^{2}}}\begin{pmatrix} {{\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta}\sin \quad \theta} - {\rho \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}\cos \quad \theta}} \\ {{{- \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta}}\cos \quad \theta} - {\rho \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}\sin \quad \theta}} \\ \rho \end{pmatrix}}}}}} & {{Formula}\quad 17} \end{matrix}$

[0074] Besides, even if the data on the corneal shape obtained through the corneal shape measurement is in the polar coordinate system form, Formula 17 cannot be directly used. As a corneal height fp (ρ,θ) is necessary to be partially differentiated with respect to θ in order to use Formula 17, a radius direction position ρ shall be constant data. This is because, since ρ is a detection edge position for the ring pattern image, it generally varies according to angles. Therefore, in order to obtain the normal vector N_(pre) of the pre-operative corneal shape, Formula 17 is used after the pre-operative cornel shape is converted to the polar coordinate system data form in which ρ does not depend on the angles by means of interpolation. Otherwise, Formula 13 is used after the pre-operative corneal shape is converted to the x-y rectangular coordinates system form.

[0075] Next, a method for calculating the shape based on the normal vector of the post-operative cornea obtained according to Formulae 5 and 7 will be described hereinafter. Assume that the normal vector N is obtained as the following components. $\begin{matrix} \begin{matrix} {\left\lbrack {{Expression}\quad 18} \right\rbrack {\overset{}{N} \equiv \begin{pmatrix} N_{x} \\ N_{y} \\ N_{z} \end{pmatrix}}} & \quad \end{matrix} & {{Formula}\quad 18} \end{matrix}$

[0076] Then, Formulae 19 and 20 hold for Formulae 17 and 18. $\begin{matrix} {\left\lbrack {{Expression}\quad 19} \right\rbrack {{{N_{x}\sin \quad \theta} - {N_{y}\cos \quad \theta}} = {{\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta}/\sqrt{\left( {\rho \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}} \right)^{2} + \left( \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta} \right)^{2} + \rho^{2}}} = {\frac{N_{z}}{\rho}\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta}}}}} & {{Formula}\quad 19} \\ {{\left\lbrack {{Expression}\quad 20} \right\rbrack  - {N_{x}\cos \quad \theta} - {N_{y}\sin \quad \theta}} = {{\rho {\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}/\sqrt{\left( {\rho \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}} \right)^{2} + \left( \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta} \right)^{2} + \rho^{2}}}} = {N_{z}\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}}}} & {{Formula}\quad 20} \end{matrix}$

[0077] Therefore, by arranging Formulae 19 and 20, formulae 21 and 22 may be obtained. $\begin{matrix} {{\left\lbrack {{Expression}\quad 21} \right\rbrack \therefore{\frac{1}{\rho}\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\theta}}} = \frac{{N_{x}\sin \quad \theta} - {N_{y}\cos \quad \theta}}{N_{z}}} & {{Formula}\quad 21} \\ {{\left\lbrack {{Expression}\quad 22} \right\rbrack \therefore\frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho}} = \frac{{{- N_{x}}\cos \quad \theta} - {N_{y}\sin \quad \theta}}{N_{z}}} & {{Formula}\quad 22} \end{matrix}$

[0078] Formula 22 is used for calculating the post-operative corneal shape. When the following recurrence formula in the meridian direction: $\begin{matrix} {\left\lbrack {{Expression}\quad 23} \right\rbrack {{{f_{p,i}\left( {\rho,\theta} \right)} = {{\left( \frac{\partial{f_{p}\left( {\rho,\theta} \right)}}{\partial\rho} \right)_{i}{\Delta\rho}} + {f_{p,{i - 1}}\left( {\rho,\theta} \right)}}},{{f_{p,0}\left( {\rho,\theta} \right)} = 0}}} & {{Formula}\quad 23} \end{matrix}$

[0079] is substituted by Formula 22 and calculated in order, the shape in the meridian direction may be obtained. Therein, Δρ is an interval of data in a radius direction perpendicular to the measurement optical axis of the measurement optical system shown in FIG. 1. Further, fp, o(ρ, θ)=0 is set so that the height on the measurement optical axis becomes 0. If the above calculation is performed for every meridian, the whole corneal shape may be constituted.

[0080] When the post-operative corneal shape is obtained as above, the ablation amount for each point may be calculated in a range on which an optical zone is set by comparing the post-operative corneal shape to the data on the corneal shape obtained from the corneal shape measurement.

[0081] Next, ablation amount calculation operation by the computer 60 will be described. When the measurement data on the objective eye refractive power and that on the corneal shape are inputted to the control unit 61, and a control program for the ablation amount calculation by the control unit 61 is executed, a simulation map display 200 as shown in FIG. 5 is displayed on the display 62. The keyboard 63 or the mouse 64 is operated to input desired refractive power, a size of a range for correction (sizes of the optical zone and a transition zone), data on the shift and tilt of an eccentric parameter to an input parameter unit 210. The calculation unit 61 a calculates the ablation amount under a condition based on the inputted data.

[0082] Displayed on the upper-left of the display 200 are a color map 201 of the pre-operative corneal shape and a color map 202 of the post-operative corneal shape. On the lower portion of the display 200, the simulated ablation amount distributions are displayed three-dimensionally. A display map 220 shows a shape map of a whole ablation amount, and displayed as height difference between the pre-operative corneal shape and the post-operative corneal shape. A display map 221 shows a shape map for which only a spherical component (axisymmetric component) is extracted. A display map 222 shows a shape map for which only a cylindrical component (symmetric component with respect to a plane) is extracted. A display map 223 shows a shape map for which only an irregular astigmatic component (irregular component) being an asymmetric component is extracted. The calculation unit 61 a obtains the irregular astigmatic component as a remaining component after the spherical and cylindrical components are subtracted from the distribution data on the whole ablation amount. When one of the display maps 220-223 is designated with the mouse, it becomes observable as a sectional shape in an arbitrary meridian direction. Further, a maximum ablation amount value is displayed under each of the display maps 220-223 as a numerical value.

[0083] Besides, assume that the shape map initially calculated and displayed is simulated without inputting the eccentric parameter, and the irregular astigmatic component displayed in the display map 223 mainly exists in an eccentric direction. As for setting of the eccentric parameter in such a case, for example, a value of the tilt is determined based on a position and a size of the eccentricity in the irregular astigmatic component so as to reduce the irregular astigmatic component. Also, a value of the shift is moved to a pupil center so as to be in alignment with an axis at the time of the ablation.

[0084] When the eccentric parameter is inputted, the calculation unit 61 a calculates the post-operative corneal shape, and a result thereof is again displayed on the display 200 as the simulation map. An operator controls the corneal ablation amount distribution with such simulation to make judgment as to whether the eccentricity of an ablation amount calculation reference axis is appropriate or not. Incidentally, the eccentricity of the reference axis for the ablation amount calculation may be automatically determined by the calculation by the calculation unit 61 a so that the asymmetric irregular astigmatic component becomes smaller than a predetermined permissible value, instead of inputting the eccentric parameter by the operator. As a result of the simulation of the ablation amount distribution, if the post-operative corneal shape is favorable, the obtained data may be recorded and stored in a floppy disk by the FDD 66 to be used by the corneal surgery apparatus.

[0085] Further, in a case where the additional ablation for correction is performed on the post-ablative eye onto which the laser beam has been irradiated with its axis deviated, it is particularly advantageous if the reference axis for the ablation amount calculation is decentered. For example, assume that, as shown in FIG. 6, the axial deviation occurs with respect to a measurement optical axis direction 150 (, which may be regarded as the approximate visual axis direction) at the time of the previous ablation, and an area 151 is ablated. If the measurement optical axis with respect to the eye is not varied from the pre-operative condition, the calculation is to be performed so as to ablate an area 152 at the time of the additional ablation. However, the actual result is that the area 151 which has already been ablated is mainly ablated, and the ablation is not allowed due to poor corneal thickness. Investigation on this problem revealed that the visual axis direction of the eye was decentered from a pre-operative direction 150 to a direction 150′ due to the axial-deviated irradiation, and the measurement was performed under such a condition. This is appreciated from the fact that the human eye is the eccentric optical system consisting of a plurality of refractive elements. In such a case, the data for decentering the reference axis for the ablation amount calculation is inputted so as to mainly ablate the area 152. The ablation on the area 152 enables correcting the previous axial-deviated irradiation, and making the shape to be approximately in the state of initially desired one. Thus, the corneal ablation amount may be reduced, and ablation applicability may be expanded even for the eye to which the ablation has not been applicable.

[0086] In the preferred embodiment as above, the apparatus for measuring the eye refractive power distribution is utilized, and the data on the corneal ablation amount is obtained based on the measurement data. However, the present invention may naturally be applied to a case where the data on the corneal ablation amount is obtained based on the measurement data on the wavefront aberration distribution. As an example of the measurement of the wavefront aberration distribution, a method may be used where an image of a light source is projected onto a fundus of the patient's eye, and a reflection light from the fundus forms an image on an image sensor via a number of microlens arrays arranged at positions conjugate with a pupil to measure a wavefront of the light refracted in a cornea of the patient's eye based on image forming information. (See U.S. Pat. No. 6,086,204 and others.) The measurement of the refractive power distribution and that of the wavefront aberration distribution are different from each other in a form of presentation of the measurement results; however, those are the same in that the optical aberration of the eye is measured. As the vector Q_(pre) (the unit direction vector of the incident light into the pre-operative cornea) expressed in Formulae 2 and 5 may be obtained via the wavefront aberration of the eye, the present invention may be applied for calculating the corneal ablation amount also in the case of the measurement of the wavefront aberration distribution.

[0087] In addition, in a method where the corneal ablation amount is calculated using an approximate expression which is known as a processing technique for an optical component, based on the measurement of the wavefront aberration of the whole eye without using the data on the corneal shape (although deviation at a peripheral part becomes greater since the approximate expression is used), the ablation amount distribution may be controlled by applying the method of decentering the reference axis for the ablation amount calculation.

INDUSTRIAL APPLICABILITY

[0088] As described above, according to the present invention, the corneal ablation amount distribution may be controlled. Therefore, unnecessary corneal ablation amount is reduced, and the ablation applicability may be expanded. 

1. A corneal-ablation-data determining apparatus for use in keratorefractive surgery for ablating a corneal surface to correct a refractive error, the apparatus comprising: input means for inputting measurement data on any one of wavefront aberration distribution and refractive power distribution of a patient's eye; decentering means for decentering a reference axis being a calculation reference for a corneal ablation amount; and calculation means for obtaining the corneal ablation amount based on the measurement data and data on eccentricity of the reference axis.
 2. The corneal-ablation-data determining apparatus according to claim 1, wherein the decentering means includes means for decentering the reference axis with respect to a measurement optical axis used for measurement of any one of the wavefront aberration distribution and the refractive power distribution, and means for inputting the data on the eccentricity.
 3. The corneal-ablation-data determining apparatus according to claim 1, wherein the input means inputs the measurement data on the refractive power distribution; and the calculation means obtains a post-operative corneal shape based on the measurement data on the refractive power distribution using a vector expression under the Snell's law so that an optical path inside the eye becomes the same before and after the surgery.
 4. The corneal-ablation-data determining apparatus according to claim 1, further comprising means for inputting data on a pre-operative corneal shape, wherein the calculation means obtains a post-operative corneal shape based on the data on the pre-operative corneal shape, the measurement data, and the data on the eccentricity of the reference axis.
 5. The corneal-ablation-data determining apparatus according to claim 1, wherein the calculation means obtains a direction vector of incident light into a post-operative cornea when the reference axis is decentered, and obtains a condition for a post-operative corneal shape as a normal vector of the cornea based on the direction vector.
 6. A corneal-ablation-data determining program for calculating a corneal ablation amount for use in keratorefractive surgery for ablating a corneal surface to correct a refractive error, the program making a computer to function as: means for inputting measurement data on any one of wavefront aberration distribution and refractive power distribution of a patient's eye; means for inputting data on eccentricity for decentering a reference axis being a calculation reference for a corneal ablation amount with respect to a measurement optical axis used for the measurement; and calculation means for obtaining the corneal ablation amount based on the measurement data and the data on the eccentricity of the reference axis.
 7. The corneal-ablation-data determining program according to claim 6, wherein the calculation means obtains a direction vector of incident light into a post-operative cornea when the reference axis is decentered, and obtains a condition for a post-operative corneal shape as a normal vector of the cornea based on the direction vector. 